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Realization of graphene moire superstructures on the surface of 4dand 5 d transition metals offers templates 
with periodically modulated electron density, which is responsible for a number of fascinating effects, 
including the formation of quantum dots and the site selective adsorption of organic molecules or metal 
clusters on graphene. Here, applying the combination of scanning probe microscopy/spectroscopy and the 
density functional theory calculations, we gain a profound insight into the electronic and topographic 
contributions to the imaging contrast of the epitaxial graphene/Ir(lll) system. We show directly that in 
STM imaging the electronic contribution is prevailing compared to the topographic one. In the force 
microscopy and spectroscopy experiments we observe a variation of the interaction strength between the tip 
and high-symmetry places within the graphene moire supercell, which determine the adsorption sites for 
molecules or metal clusters on graphene/Ir(lll). 
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Graphene layers on metal surfaces have been attracting the attention of scientists since several decades, 
starting from middle of the 60s, when the catalytic properties of the close-packed surfaces of transition 
metals were in the focus of the surface science research 1 " 4 . The demonstration of the fascinating electronic 
properties of the free-standing graphene 5 ' 6 renewed the interest in the graphene/metal systems, which are 
considered as the main and the most perspective way for the large-scale preparation of high-quality graphene 
layers with controllable properties 7 " 10 . For this purpose single-crystalline as well as polycrystalline substrates of 
3d - 5d metals can be used. 

One of the particularly exciting questions concerning the graphene/metal interface is the origin of the 
bonding mechanism in such systems 2 " 411 . This graphene- metal puzzle is valid for both cases: graphene 
adsorption on metallic surfaces as well as for the opposite situation of the metal deposition on the free- 
standing or substrate-supported graphene. In the later case the close-packed surfaces of 4d and 5d metals are 
often used as substrates 12 " 14 . A graphene layer prepared on such surfaces, i. e. Ru(OOOT) 15 " 17 , Rh(lll) 1418,19 , 
Ir(lll) 20 " 24 , or Pt(lll) 25 ' 26 , forms so-called moire structures due to the relatively large lattice mismatch 
between graphene and metal substrates. As a consequence of the lattice mismatch the interaction strength 
between graphene and the metallic substrate is spatially modulated leading to the spatially periodic electronic 
structure. Such lateral graphene superlattices are known to exhibit selective adsorption for organic mole- 
cules 27 or metal clusters 28 . Especially, the adsorption of different metals - Ir, Ru, Au, or Pt - on graphene/ 
Ir(lll) has been intensively studied showing a preferential nucleation around the so-called FCC or HCP 
high-symmetry positions within the moire unit cell 13 ' 29 . In the subsequent works 29 ' 30 this site -selective adsorp- 
tion was explained via local sp 2 to sp 3 rehybridization of carbon atoms with the bond formation between 
graphene and the cluster. However, a fully consistent description of the local electronic structure of gra- 
phene/Ir(lll), the observed imaging contrast in scanning probe experiments and the bonding mechanism of 
molecules or clusters on it is still lacking, motivating the present research. 

Here we present the systematic studies of the graphene/Ir(lll) system by means of density functional theory 
(DFT) calculations and scanning tunnelling and atomic force microscopy (STM and AFM) performed in con- 
stant current/constant frequency shift (CC/CFS) and constant height (CH) modes. The obtained results for the 
graphene/Ir( 111) system allow to separate the topographic and electronic contributions in the imaging contrast in 
STM and AFM as well as to shed light on the spatially modulated interaction between graphene/Ir(lll) and the 
metallic STM/AFM tip, which is of paramount importance for the understanding of absorption of metals on top 
of graphene/Ir(lll) as well as of similar graphene-metal systems. 
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Results 

The unit cell of graphene on Ir(lll) is shown in Fig. 1(a) with the 
corresponding high- symmetry local arrangements of carbon atoms 
above Ir layers marked in the figure: ATOP (circles), FCC (squares), 
HCP (down-triangles), and BRIDGE (stars). The DFT-D2 optimized 
local distances between graphene and Ir(l 1 1) are 3.27 A (HCP), 3.28 
A (FCC), 3.315 A (BRIDGE), 3.58 A (ATOP). This result is very close 
to the recently published equilibrium structure for this system 21 . The 
similar distances for HCP, FCC, and BRIDGE positions can be related 
to the fact that in all these cases one of the carbon atoms in the 
graphene unit cell is placed above Ir(S) atom defining the local inter- 
action strength for the particular high-symmetry position. The 
obtained distances are very close to those between carbon layers in 
pure graphite and this was explained by a binding interaction domi- 
nated by van der Waals effects between graphene and Ir(l 1 1) that is 
modulated by weak bonding interactions at the FCC and HCP places 
and anti-bonding chemical interaction around ATOP positions 21 . As 
a result a small charge transfer from the graphene n states on Ir 
empty valence band states is detected in calculations and the Dirac 
point is shifted by ~ 100 meV above the Fermi level (E F ) that is close 
to the data obtained by photoelectron spectroscopy 31 [see Fig. 2(c) 
and discussion below]. 

The graphene/Ir(l 1 1) system is a nice example of the moire struc- 
ture 32 " 37 , which is easily recognisable in LEED and on the large scale 
STM images shown in Fig. 1(b). The extracted lattice parameter of 
this structure from LEED and STM is 25.5 A and 25.2 A, respectively, 
that is in good agreement with previously published data 32 . In most 
cases, for the typical bias voltages used in STM imaging, the gra- 
phene/Ir(lll) structure is imaged in the so-called inverted contrast 
[Fig. 1 (c,d)] 32 , when topographically highest ATOP places are imaged 
as dark and topographically lowest FCC and HCP ones as bright 
regions. This assignment was initially done in Ref. [32], where the 
STM topography of neighbouring regions of the clean and graphene- 
covered Ir ( 111) surface were imaged. In our studies we perform the 
"on-the-fly" switching between CC STM and CFS AFM imaging 
during scanning [Fig. 1(d)] where we observe the inversion of the 
topographic contrast z(x, y) (see also Fig. SI of the supplementary 
material). This becomes clearly evident around areas marked by 
arrows in Fig. 1(d) where the darkest contrast in CC STM becomes 
the brightest one in CFS AFM for the ATOP position. In the later case 
the darkest areas correspond to HCP sites that correlates with the 
calculated height changes between lowest and highest carbon posi- 
tions. The pronounced difference between the FCC and HCP areas 
can be explained by fact that additional bias voltage was applied in 
order to increase the atomically-resolved contrast. The present result 
is opposite to the observations published earlier in Ref. [37] and 
where no contrast inversion between STM and AFM was observed. 
Our further analysis demonstrates that the imaging contrast strongly 
depends on the experimental conditions (tunnelling bias and cur- 
rent, frequency shift, distance from the sample), but performed 
switch "on the fly" between STM and AFM demonstrates the pos- 
sibility to directly trace the changes between STM and AFM pictures. 

The inversion of the imaging contrast was also detected in CC 
STM images when the bias voltage is changed from —0.5 V to 
— 1.8 V during scanning [Fig. 2(a)]. The images taken at positive 
voltages are similar to the lower part of the figure Fig. 2(a) (see 
Fig. S2 of the supplementary material). The simulations of STM 
images using the experimental conditions show very good agree- 
ment with the obtained data [Fig. 2(a,b) and Fig. S2]. The pre- 
sented STM results demonstrate the big difference in the local 
electronic structure for high-symmetry positions of carbon atoms 
on Ir(lll). That can point out on the difference in the local 
adsorption strength for these places. 

In order to prove this assertion we have performed the force 
microscopy and spectroscopy experiments on graphene/Ir(lll). 
Fig. 3 shows (a,b) the frequency shift of an oscillating scanning sensor 




• C (graphene) • Ir (S) O Ir (S-1) O Ir (S-2) 




Figure 1 | (a) Crystallographic structure of (10 X 10) graphene on (9 X 9) 
Ir( 111). The high-symmetry places are marked by circle, rectangle, down- 
triangle, and stars for ATOP, FCC, HCP, and BRIDGE positions, (b) Large 
scale (162 X 55 nm 2 ) STM image of graphene/Ir( 111). The inset shows the 
corresponding LEED image obtained at 71 eV. (c) Atomically resolved 
STM image of graphene/Ir(lll) showing the inverted contrast, (d) A 
combined STM/ AFM image with a switching "on-the-fly" between CC 
STM and CFS AFM imaging during scanning. The scan sizes are 5.2 X 
5.2 nm 2 (U T = +550 meV, I T = 540 pA) and 10 X 10 nm 2 (STM: U T = 
+461 meV, I T = 7.6 nA, AFM: Af= -675 mHz) for (c) and (d), 
respectively. 

as a function of a distance from the sample, Af(d), and (c) the cor- 
responding tunnelling current, 1(d), measured in the unit cell of 
graphene/Ir(lll) along the path marked in the CC STM image 
shown as an inset of (a). During these measurements the tunnelling 
current was used for the stabilisation of the feedback loop, allowing to 
determine the relative z-position of the A/and I curves. Following the 
sequence: F z (d) = —dE(d)ldd, Af(d) = —f 0 /2k 0 'dF z (d)/dd, where 
E(d) and F z (d) are the interaction energy and the vertical force 
between a tip and the sample, respectively, / 0 and k 0 are the resonance 
frequency and the spring constant of the sensor, correspondingly, we 
can separate repulsive, attractive, and long-range electrostatic and 
van der Waals contributions in AFM imaging. The imaging contrast 
at distances more than 5 A from the surface is defined by long-range 
van der Waals interactions which are insensitive to the atomically 
resolved structure of the scanning tip and the sample. The short- 
range chemical interactions around the minimum of the A/curve are 
dominated either by repulsive (left-hand side) or attractive (right- 
hand side) forces and give the atomically-resolved site- selective 
chemical contrast. 
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Figure 2 | Experimental (It = 15 nA) (a) and calculated (b) CC STM images of graphene/Ir( 111) obtained at —0.5 V (bottom) and —1.8 V (top) of the 
bias voltage, (c) Carbon-site projected density of states calculated for all high-symmetry places of graphene/Ir(l 1 1) (zoom around E F is shown as an 
upper inset). Difference electron density, Ap(r) = p gr /ir(m)( r ) ~ PiM ~ P g M> plotted in units of e/A 3 is shown as lower inset. 



The presented Af(d) curves [Fig. 3(a,b)] clearly show the site- 
selective interaction in the unit cell of the graphene/Ir(lll) system. 
First of all, the absolute value of the maximal frequency shift is higher 
for the FCC and HCP areas compared to ATOP ones indicating the 
stronger interaction of graphene with the tungsten STM/AFM tip for 
the former carbon positions. We would like to note, that the same 
trend will be observed for the adsorption of metal atoms or clusters 
on graphene/Ir(lll). Secondly, the difference in the positions of 
A/(d)-curve minima is 0.097 nm [Fig. 3(b)] compared to the value 
of graphene corrugation of 0.067 nm obtained from the CC STM 
image [inset of Fig. 3(a)], also pointing on the difference in the 
interaction strength for two carbon positions, ATOP and FCC. 

The measured Af(d) curves were used to track and explain the 
contrast obtained in CH AFM experiments (Fig. 4 and Fig. S4 in 
supplementary material). In this experiments the z- coordinate of 
the oscillating sensor was fixed during scanning and Af(x, y) and 
I(x, y) maps were collected. As can be clearly seen from the presented 
results the imaging contrast in CH Af images is changing as a func- 
tion of z- coordinate and it is fully inverted for two limit z-positions 
[di = —0.033 nm and d 2 = +0.227 nm, see Fig. 3(b)] used for the 
imaging (upper row in Fig. 4). In the region of distances between the 
tip and the graphene/Ir(lll) sample, where chemical forces deter- 
mine the interaction, on the right-hand side of the Af curve the 
interaction is attractive giving the more negative frequency shift 
for the ATOP positions as located closer to the oscillating tip com- 
pared to the FCC and HCP areas. On the other side, on the left-hand 



side of the curve the repulsive interaction starts to contribute in the 
interaction giving the smaller frequency shift for ATOP compared to 
other positions that can again be explained by the different distances 
between tip and the sample for different high-symmetry positions as 
well as by the different interaction strength at these places. At the 
same time the contrast for the current map I(x,y) is always the same, 
inverted contrast (lower row in Fig. 4). Taking into account that 
during CH imaging at + 50 meV in Fig. 4 the imaging contrast for 
I(x, y) is inverted due to the higher DOS for the FCC and HCP 
positions and that variation of the distance for CH imaging and 
CC STM imaging (when U T is changed from —0.5 V to —1.8 V) is 
nearly the same (see next section for discussion of STM simulations), 
one can separate the topographic and electronic contributions into 
imaging at different biases and distances. 

Here we would also like to note that the change of the bias voltage 
during CH AFM imaging does not lead to any changes in the imaging 
contrast for Af (Fig. 5). The first row shows the two small -scanning 
range atomically resolved Af(x, y) images acquired on the same place 
of the graphene/Ir(l 1 1) sample with opposite signs for the bias volt- 
age: the imaging contrast is the same with the slight variation of the 
imaging scale that can be explained by the small drift of the oscil- 
lating tip. The first look on the I(x,y) map (lower row) might give an 
impression that the contrast is fully inverted. However, the absolute 
value of the tunnelling current is the same and only more negative 
values of the tunnelling current are shown as darker areas in the 
image for current. 
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Figure 3 | Frequency shift (a,b) of the oscillating KolibriSensor and the corresponding tunnelling current (c) as a function of the relative distance d 
between tip and the graphene/Ir(l 1 1) sample. These curves were acquired with respect to the set-point for the CC STM imaging. Inset of (a) shows the 
corresponding STM image (U T = +50 mV, I T = 400 pA) with the path where Af and I data were measured. 



Discussion 

The graphene/Ir(lll) system was studied with DFT methods in 
order to deeply understand the effects of the electronic structure 
on the adsorption properties of this system. Here we present the 
comparison between experimental and theoretical STM data. The 
STM images are calculated using the Tersoff-Hamann formalism 38 , 
in its most basic formulation, approximating the STM tip by an 
infinitely small point source 19 ' 39,40 . In these simulations the constant 
current condition was fulfilled that leads to the increasing of the 
distance between the sample and a tip from 2.50 A for U T = 
-0.5 V to 3.21 A for U T = -1.8 V [Fig. 2(b)]. Taking into account 



that at negative bias voltages the electrons tunnel from occupied 
states of the sample into unoccupied states of the tip and that tun- 
nelling current is proportional to the surface local density of states 
(LDOS) at the position of the tip we can clearly identify the states in 
the valence band of graphene/Ir(lll) which are responsible for the 
formation of the imaging contrast (see Fig. S3 of the supplementary 
material for the identification of the local valence band states states). 
[The surface Brillouin zone (BZ) of the graphene/Ir(lll) system is 
approximately 10 times smaller compared to BZ of free-standing 
graphene giving the maximal wave-vector of electrons of k\\ ~ 0.17 
A" 1 for the iC-point, i. e. in this case one has to consider a small region 



z=d u U T =50mV z=d 2 , U T =50mV 

^_^+1.02Hz -i 




Figure 4 | Constant height images of graphene/Ir(l 1 1), A/(x, y) (upper row) and I(x, y) (lower row), obtained at two different heights d x and d 2 of the 
KolibriSensor above the surface (see Fig. 3). The size of the scanning area is 10.5 X 10.5 nm 2 . 



SCIENTIFIC REPORTS | 3 : 1072 | DOI: 1 0.1 038/srep01 072 



4 



U T = +400 mV U T = -400 mV 

-1.00Hz -1.19Hz 



2.62 Hz 




238 pA 



277 pA 



385 pA 



Figure 5 | Atomically- resolved constant height images of graphene/Ir(lll), A^x, y) (upper row) and I(x, y) (lower row), collected at opposite signs of 
bias voltages. The size of the scanning area is 4.3 X 4.3 nm 2 . 



in the reciprocal space around the T-point. In this case we can assign 
the features in the tunnelling current as originating from the corres- 
ponding peaks in DOS of the whole system (in general, the tunnelling 
current depends on the k\\ of electronic states participating in the 
tunnelling process and the electronic states with a large parallel com- 
ponent of the wave vector give a small contribution to the current)]. 
Therefore, the so-called inverted contrast for the small bias voltages, 
negative or positive, can be explained by two features in LDOS for the 
carbon atoms in the FCC position located at E — E F = — 0.53 eV and 
E — E F = +0.32 eV, respectively [Fig. 2(c)] . The next two maxima in 
LDOS for C ATOP positions located at E - E F = -0.95 eV and E - 
E F = — 1.49 eV [Fig. 2(c)] explain the fact that these high-symmetry 
places become brighter for higher bias voltages in experimental and 
simulated CC STM images [Fig. 2(a,b)]. 

The calculated LDOS for graphene/Ir(lll) compared with DOS 
for the free-standing graphene [(10 X 10) unit cell] are shown in 
Fig. 2(c) together with the electron density difference distribution 
presented as lower inset. The obtained DOS results demonstrate the 
slight difference in the position of main features (71-states at E — E F ~ 
6.1 eV, cr-states at E — E F ~ 3.2 eV, and the M-point- derived Kohn 
anomaly at E — E F ~ 2.2 eV) as well as the Dirac cone between two 
systems, that is assigned to thep-doping of graphene on Ir(l 1 1). The 
DOS features responsible for the formation of the STM contrast can 
be considered as a result of the hybridisation between graphene n 
states and the valence band states of Ir around the particular high- 
symmetry adsorption positions (see Fig. S3 of the supplementary 
material for the identification of the local valence band states states). 
The trace of this hybrid states can be identified in the full photoemis- 
sion map of the graphene/Ir(lll) system presented in Ref. [11]: the 
calculated DOS peaks can be identified as weak "flat" hybrid bands in 
experimental band structure along the K — M direction. However, as 
was discussed earlier the folding of the bands leads to the fact that 
these bands will appear around the Y point and can effectively con- 
tribute to the tunnelling and can be responsible for the contrast 
inversion observed in the experiment. Here we would like to note, 
that one can also consider the increased LDOS around E F for FCC, 
HCP, and BRIDGE positions as a hint indicating that these places can 
be nucleation centers for adsorbed atoms 1213 . 

In order to get better insight in the results obtained during our 
AFM experiments (Figs. 3, 4, 5) we performed simulations of these 



data in the framework of DFT formalism (we would like to note that a 
qualitative description of the buckled graphene systems was per- 
formed in Ref. [41], where the model Lennard- Jones potential was 
used to model tip-sample interaction). Due to the large size of the 
unit cell of the graphene/Ir(lll) system we employed two 
approaches. In the first one, which requires less computational 
resources, the tip-sample force is expressed as a function of the 
potential V fs (r) on the tip due to the sample: F ts cc— V[|VV7 5 
(r)| 2 ] 42 . However, this method does not take into account the geo- 
metrical and electronic structure of the scanning tip and thereby does 
not allow to get absolute values for the force or frequency shift and 
the correct distance between the tip and the sample, and only qual- 
itative result for the attractive region of the interaction can be 
obtained. The result of such simulations for graphene/Ir(lll) is 
shown in Fig. 6(a), which is in rather good agreement with experi- 
mental data: the regions around ATOP positions are imaged as dark 
compared to other high-symmetry cites imaged as brighter contrast. 
Unfortunately, the information about repulsive region of the forces 
cannot be obtained from such calculations. However, if the model 
corrugation -dependent repulsive potential is added, then the result- 
ing simulated Af curve reproduces qualitatively the experimentally 
obtained results. 

For the purpose to reproduce our data in a more quantitative way, 
the interaction between the tip and the graphene/Ir(l 1 1) system was 
simulated within the second approach, where W-tip is approximated 
by the 5 -atom pyramid as shown in Fig. 6(b). The results of these 
calculations are compiled in Fig. 6(c-e). The interaction energy (sys- 
tem was rigid without allowing to relax) between model W-tip and 
the surface was calculated for two limiting places of graphene/ 
Ir(lll), ATOP and FCC. The calculated points are shown by filled 
rectangles and circles in Fig. 6(c) for the FCC and ATOP positions, 
respectively. The Morse potential was used to fit the calculated data 
as the most suitable for the graphene-metal systems 43,44 . The resulting 
curves are shown by solid lines in the same figure. The force and the 
frequency shift curves were calculated on the basis of the obtained 
interaction energy curves according to formulas presented earlier. 
The results are shown in Fig. 6(d) and (e), respectively. 

The interaction energy curves [Fig. 6(c)] are very similar for the 
two high- symmetry areas and shifted with respect to each other by 
0.40 A reflecting the height difference between ATOP and FCC 
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Figure 6 | (a) Simulated CH AFM image of the graphene/Ir( 111) system according to the approach suggested in Ref. [42] . (b) Schematic representation of 
the geometry of the 5-atom-W-tip/graphene/Ir(lll) system used in simulation of the frequency shift curves. Interaction energy (c) and the force (d) 
between 5-atom W-tip and graphene/Ir(l 1 1) calculated for two limit places ATOP and FCC. (e) Frequency shift as a function of distance between the 
model 5-atom-W-tip and graphene/Ir(lll). 



regions. The maximal absolute value of energy is slightly higher for 
FCC by 0.015 eV. However, our calculations do not take into account 
the relaxation of the system during the tip -sample interaction. 
Therefore the actual values can be different. Calculation of the force 
and the frequency shift via derivation of the respective curves leads to 
the clear discrimination in the value of the extrema as well as to the 
further separation of their positions and for the A/curves the extrema 
are separated by 0.45 A. The resulting values for the frequency shift 
for our model W-tip are approximately 2 times larger compared to 
experimental values. Nevertheless, the shape and the trend for the 
inversion of the imaging contrast are clearly reproducible in our 
calculated data. The following reasons can explain the difference 
between experiment and theory: (i) small W- cluster modelling the 
actual W-tip, (ii) system was not relaxed during calculations, (hi) the 
actual atomic structure of the tip during experiment is unknown, (iv) 
small uncertainty in the tip position during measurements of Af 
curves. However, in spite of these simplifications 45 , the presented 
calculation results allow us to confirm and better understand the 
data obtained during the AFM measurements. 

In conclusion, we performed the systematic investigations of the 
geometry, electronic structure and their effects on the observed 
imaging contrast during STM and AFM experiments on graphene/ 
Ir(lll). Our results obtained via combination of DFT calculations 
and scanning probe microscopy imaging allow to discriminate the 
topographic and electronic contributions in these measurements and 
to explain the observed contrast features. We found that in STM 
imaging the electronic contribution is prevailing compared to the 
topographic one and the inversion of the contrast can be assigned to 
the particular features in the electronic structure of graphene on 
Ir(lll). Contrast changes observed in constant height AFM mea- 
surements are analyzed on the basis of the energy, force, and fre- 
quency shift curves reflecting the interaction of the W-tip with the 



surface and are attributed to the difference in the height and the 
different interaction strength for high-symmetry cites within the 
moire unit cell of graphene on Ir(lll). The presented findings are 
of general importance for the understanding of the properties of the 
lattice-mismatched graphene/metal systems especially with regard to 
possible applications as templates for molecules or clusters. 

Note added after submission. During consideration of the present 
manuscript similar AFM results on graphene/Ir(lll) were 
published in Ref. [46] . 

Methods 

DFT calculations. The crystallographic model of graphene/Ir(lll) presented in 
Fig. 1(a) was used in DFT calculations, which were carried out using the projector 
augmented plane wave method 47 , a plane wave basis set with a maximum kinetic 
energy of 400 eV and the PBE exchange-correlation potential 48 , as implemented in 
the VASP program 49 . The long-range van der Waals interactions were accounted for 
by means of a semiempirical DFT-D2 approach proposed by Grimme 50 . The studied 
system is modelled using supercell, which has a (9 X 9) lateral periodicity and 
contains one layer of (10 X 10) graphene on a four-layer slab of metal atoms. Metallic 
slab replicas are separated by ca. 20 A in the surface normal direction. To avoid 
interactions between periodic images of the slab, a dipole correction is applied 51 . The 
surface Brillouin zone is sampled with a (3 X 3 X 1) A;-point mesh centered the T 
point. 

STM and AFM experiments. The STM and AFM measurements were performed in 
two different modes: constant current (CC) or constant frequency shift (CFS) and 
constant height (CH). In the first case the topography of sample, z(x, y), is studied 
with the corresponding signal, tunnelling current (Jj) or frequency shift (Af), used as 
an input for the feedback loop. In the latter case the z-coordinate of the scanning tip is 
fixed with a feedback loop switched off that leads to the variation of the distance d 
between tip and sample. In such experiments Ij{x, y) and Af(x, y) are measured for 
different z-coordinates. The STM/ AFM images were collected with SPM Aarhus 150 
equipped with KolibriSensor™ from SPECS 19 ' 52 with Nanonis Control system. In all 
measurements the sharp W-tip was used which was cleaned in situ via Ar + - 
sputtering. In presented STM images the tunnelling bias voltage, U T , is referenced to 
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the sample and the tunnelling current, I T , is collected by the tip, which is virtually 
grounded. During the AFM measurements the sensor was oscillating with the 
resonance frequency of/ 0 — 999161 Hz and the quality factor of Q = 45249. The 
oscillation amplitude was set to A = 100 pm or A = 300 pm. 

Preparation of graphene/Ir(lll). The graphene/Ir(lll) system was prepared in 
ultrahigh vacuum station for STM/AFM studies according to the recipe described in 
details in Ref. [32] via cracking of ethylene: T = 1100°C, p = 5 X 10" 8 mbar, t = 
15 min. This procedure leads to the single-domain graphene layer on Ir(l 1 1) of very 
high quality that was verified by means of low-energy electron diffraction (LEED) and 
STM. The base vacuum was better than 8 X 10" 11 mbar during all experiments. All 
measurements were performed at room temperature. 
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